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Abstract 

We point out that single electron charging effects such as Coulomb Blockade (CB) and high- 
bias staircases play a crucial role in transport through single ultrashort molecules. A treatment of 
Coulomb Blockade through a prototypical molecule, benzene, is developed using a master-equation in 
its complete many-electron Fock space, evaluated through exact diagonalization or full Configuration 
Interaction (CI). This approach can explain a whole class of non-trivial experimental features including 
vanishing zero bias conductances, sharp current onsets followed by ohmic current rises, and gateable 
current levels and conductance structures, most of which cannot be captured even qualitatively within 
the traditional Self Consistent Field (SCF) approach coupled with perturbative transport theories. By 
comparing the two approaches, namely SCF and CB, in the limit of weak coupling to the electrode, 
we establish that the inclusion of strong-correlations within the molecule becomes critical in addressing 
the above experiments. Our approach includes on-bridge-correlations fully, and is therefore well-suited 
for describing transport through short molecules in the limit of weak coupling to electrodes. 

Index Terms 
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Molecular conduction, Self Consistent Field, Non Equilibrium Green's Function, Coulomb Block- 
ade, Configuration Interaction. 

I. Introduction 

Theoretical calculations on single molecule conduction have typically employed coherent Non- 
Equilibrium Green's function (NEGF) theories ("Landauer limit") [1], [2] coupled with Self 
Consistent Fields (SCF) to describe charging effects. Though fairly successful in describing many 
aspects of single molecule conduction [3]-[7], there have been important discrepancies between 
theory and experiment [8]. The most common ones include poor match between theoretical and 
experimental current levels and zero-bias currents [3], [4], [8]. It was also pointed out in [9] that 
a whole class of experimental I-V's show features, which cannot be captured even qualitatively 
using an SCF theory. Charging energies of short molecules ( 3 eV for benzene) are often larger 
than their electrode coupling (< 0.2 eV for benzene di-thiol on gold), and thus could be in the 
Coulomb Blockade (CB) regime where single electron charging effects could dominate. It is thus 
debatable whether it is better described as a quantum wire in the SCF regime, or as a quantum dot 
array in the Coulomb Blockade (CB) regime. Nevertheless the wisdom of SCF approaches must 
be scrutinized especially for conduction through shorter molecules. The purpose of this paper is 
to present a Coulomb Blockade approach to molecular conduction using a benzene molecule as 
prototype, and establish it as a different viewpoint from the conventional NEGF-SCF treatment. 
Furthermore features obtained via the CB approach can semi-quantitatively explain several non- 
trivial features commonly observed [10]— [14] in experiments. 

It is common to distinguish between two regimes of transport: a) an SCF regime where the 
dominant energy scale is the contact coupling, allowing for fractional charge transfer through 
the system; and b) a Coulomb Blockade (CB) regime where the dominant energy scale is the 
single electron charging, leading to integral charge transfer. In the SCF regime the description 
of transport via non-interacting single particle energy levels can be justified. In this limit, it 
is common to use the SCF-NEGF scheme that takes charging effects into account. Here the 
molecular Hamiltonian is described by a set of single particle levels, which are coupled to 
reservoirs through their self energies. The electron interactions are taken into account using SCF 
schemes as shown in the block diagram in Fig. la. All quantities in the NEGF formalism are 
matrices of dimension N x N, N being the number of single particle basis functions used. This 
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Fig. 1. Our system is a benzene molecule coupled to metallic contacts. Single molecule transport calculations typically employ 
the NEGF-SCF prescription. The Block diagram depicts the basic scheme. While quantities such as the Hamiltonian are in the 
one-electron space of dimensions N x N, N being the number of basis functions. Our Coulomb Blockade description involves 
the full many electron Fock space of dimensions 2^ x 2 N as shown in b) using a single spin-degenerate level as an example. 
The use of full many-electron space captures the correlations exactly within the framework of the given N x N one-electron 
Hamiltonian. 



allows for an accurate description of quantum chemistry of both the isolated molecule and its 
bonding to the contacts [15]. In the CB regime with weak contact coupling, charging effects 
dominate, and the use of single particle basis sets may be questionable. In such cases, it may 
be preferential to employ a multi-electron description or Configuration Interaction (CI) where 
feasible. The central quantities in this CI method are now matrices of dimension 2^ x 2 N , 
thereby accounting for strong interaction accurately. The weakly coupled contacts are treated 
perturbatively using transition rates between states differing by a single electron [5]. It is 
interesting to note that most theoretical efforts in molecular conduction have been in the SCF 



February 6, 2008 



DRAFT 



4 



regime, while energy scales favor the CB regime. Our paper is thus a concrete attempt towards 
CI based transport. 

This paper is organized as follows: we begin by defining an appropriate many -body Hamil- 
tonian for Benzene whose parameters are benchmarked based on well-established mean-field 
techniques. We then illustrate how a CB treatment is conceptually different from the standard 
SCF treatment in the weak coupling limit, not only under non-equilibrium conditions, but even 
under equilibrium conditions. We then point out the importance of inclusion of excited states in 
transport, that naturally arise within our CI approach. The progressive access of these excited 
states leads to transport signatures under various non-equilibrium conditions. Before we conclude, 
a few CB fits to experimental data are presented in support of our analysis. 

II. The Model Hamiltonian and Equilibrium Properties 

An appropriate model Hamiltonian is usually described with an adequate basis set. In this 
paper, we use a tight binding Hamiltonian with one p z orbital per site to describe our CI based 
scheme. Although this generates just a minimal 6x6 single particle basis set, its many-electron 
space is 2 12 x 2 12 in size. Besides, our objective here is to describe the CI approach for transport 
and compare it with the SCF approach for the same Hamiltonian. Better quantum chemical 
descriptions within the CI approach can be achieved by starting with a reduced but more accurate 
one-particle Hamiltonian, but we leave these for future work. 

One begins with the model Hamiltonian in second quantized notation: 

+ U aa n aa n a „ + - Y U a pn a n p , (1) 

where a, (3 correspond to the orbital indices of the frozen p z orbitals for carbon sites on the 
Benzene ring,and a, a represent a particular spin and its reverse. In connection to its equilibrium 
configuration, it is more convenient to work with onsite energies e defined as: 

e a = e« + U aa {n aS ) + - ^ t^W), (2) 

where e a 's denote the mean-field on-site energies in the equilibrium charge neutral configuration 
of the molecule and (n) represents its mean-field value. Now the model Hamiltonian is simply 
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re- written as: 

H = E ^n a + E tapcicp 

+ 9 E - M)(n p - (n^)). (3) 

The mean-field Hamiltonian derived from the above Hamiltonian, is: 

h = J2 + E ^/3 c l c /3 + ^fa ^ ( 4 ) 

where 

f/f Q CF = U aa (n a " + ^ E ^(n, " (*/>»■ (5) 

is the Self Consistent Field, the calculation of (n a ) performed self consistently with the one 
electron Hamiltonian h. In the following sections, we derive appropriate parameters e, t and U 
for benzene, to describe the two different approaches i.e., the CI (Eq. O and the SCF approaches 
(Eq. 13), and compare them in parallel in the case of both equilibrium and non-equilibrium 
conditions. 

Fig. 2(a) shows the selection of mean field on-site energies e a and hopping parameter t a p, by 
comparing the eigen-energies of our model SCF Hamiltonian (Eq. |4j) with the frontier orbitals 
within the local density approximation (LDA) in the 6-3 lg basis set, shown in the left and right 
sections of Fig. 2a respectively. The carbon-carbon hopping term t a p = —2.0eV has been used 
from already tabulated data [16], which yields e a = —AA2eV for the above fit. Note that the 
Highest Occupied Molecular Orbital (HOMO) levels and Lowest Unoccupied Molecular Orbital 
(LUMO) levels are doubly degenerate in our model tight binding Hamiltonian as well as in the 
LDA basis set. 

A. Equilibrium Electron Number v/s chemical potential: Choosing Charging parameters 

A distinguishing aspect of CB is the abrupt charge addition as opposed to a gradual one in an 
SCF calculation shown in Fig. 2(b). This fact is readily seen in the figure, in which the SCF and 
CB calculations are presented using the one electron and many-electron Hamiltonians described 
by Eq. |3] and |4] respectively. In the weak coupling limit, the CB result is more physical, and 



February 6, 2008 



DRAFT 



6 




Fig. 2. Model Hamiltonian and Equilibrium Properties, (a) Selection of on-site e a and hopping parameter t a p. Comparison 
of our model Hamiltonian levels with frontier LDA/6-31g levels. Parameters are fixed based on a close match between the 
doubly degenerate HOMO and LUMO levels and singly degenerate HOMO-1, LUMO+1 levels, (b) Charging parameter matched 
according to a consistent Restricted SCF based N — fj, plot (shown continuous line). Total energy based Many-Body calculation 
(shown dotted line) as well as RSCF calculation is consistent with Gaussian based calculation [22]. c) One particle spectral 
function shows peaks at the energy levels of the single particle Hamiltonian h. d) Lehmann spectral function evaluated via many- 
electron spectrum yields many more spectral peaks corresponging to removal (addition) of electrons from the neutral ground 
state into various charge configurations (excitations) of singly charged species. Notice that the IP-EA and HOMO-LUMO gaps 
are equal to the corresponding charge-stability plateaus N = No for many-body and SCF calculations shown in b). 
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the SCF calculation does not do justice to this integer charge transfer. However, schemes such 
as self interaction correction could be introduced within the one-electron Hamiltonian [17]— [19] 
to incorporate this. But it turns out that even such schemes may not capture non-equilibrium 
correctly [9]. It is however expected that as coupling strength to the electrodes is increased, the 
electron transfer resembles the SCF result. There is as yet no clear formalism that addresses [20] 
this crossover, even in the equilibrium case although the two opposite limits namely SCF and CB 
are well understood. While the two limits can individually be handled by perturbative expansion 
in the small parameters U/Y and Y/U, U being the single electron charging energy and Y being 
the level broadening, the intermediate regime is hard to handle owing to the non-existence of a 
suitable small parameter or 'fine structure constant' for transport. 

The plateau in the charge addition diagram N versus fx, in which the electron number is 
stabilized, spans the HOMO-LUMO gap in the SCF case and the Ionization potential-Electron 
Affinity (IP-EA gap) in the CB calculation. The IP (EA) is defined as the energy when an electron 
can be removed (added) to the neutral molecule carrying N electrons. This occurs when the 
chemical potential fx equals the energy difference between ground states differing by an electron 
number fx = Eq° -E^ " 1 for IP, and fx = E^ 0+1 -E^° for EA. The situation is however different 
in the case of SCF. Here the charge transfer dictated by a self consistent potential (Eq. |5J is 
gradual, in which two electrons are transferred adiabatically over a span of 2U corresponding to 
the removal of two electrons. This is usually referred to as the restricted SCF (RSCF mentioned 
in Fig. 2b). Most SCF calculations in the literature [21] employ different variants of this scheme. 
There are also spin unrestricted SCF techniques [17]— [19] which take into account the abrupt 
charge transfer in a weakly coupled system, due to self-interaction correction, but it is not yet 
clear whether they work out of equilibrium [9], [23]. 

One expects that the IP occurs roughly midway during the gradual charge removal in the 
RSCF scheme [20], [21]. We use this fact to estimate our charging parameters U a p, with the aid 
of a Gaussian-98 based calculation for the equilibrium electron number v/s chemical potential 
(N — fx), published elsewhere [22]. The calculation corresponding to the equilibrium number of 
No = 42 maps onto our model calculations for No = 6, focussing thus on the frontier orbitals 
and ignoring the inner core that is frozen in our estimate for e a . By implementing a Restricted 
SCF scheme using Eq. 13 within in our model Hamiltonian, we obtain a close match of the N — fx 
plots in the range between N and N — 1 in comparison with the Gaussian-98 calculation in [22]. 
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Using an estimate of the onsite charging U aa , we calculate U a p using the Matago-Nishimoto 
approximation: 



e 2 



U a /3 — ; ^2 > (6) 

47re o r ^ + U^+U^ 

where r a p is the inter carbon distance in benzene. In each case, evaluation of n a is done self- 
consistently using an equilibrium value N = 6, and (n aa ) = \. Using exact eigen-energies of 
the many-electron Hamiltonian Eq.|3]wifh the above parameters, an N — p, calculation using these 
total energies (shown dotted red in Fig. 2b) is in excellent agreement with respect to Gaussian 
calculations in [22]. Note that the fi = IP in Fig. 2b occurs midway between N = N and 
N = N — 2 in the RSCF charging diagram. It is worth mentioning that the many -body calculation 
presented in this figure takes all correlation energies into account and is the exact ground state 
energy within our defined model Hamiltonian. 



B. Equilibrium Spectral function 

Conduction through molecules via molecular orbitals is well understood in the SCF picture 
[24]. In the strongly coupled regime (most appropriate for an SCF treatment), fractional charge 
transfer occurs, and Density of States (DOS) is evaluated at equilibrium [1] in order to capture the 
effect of the strong coupling with contact. An interplay of molecular DOS and charging treated 
self-consistently determines the non-equilibrium response (current-voltage or I-V characteristics). 
The density of states calculated from the one-electron Green's function [1] in Fig. 2c shows 
peaks at the single electron eigen spectrum. As the coupling to electrodes gets stronger, the 
single electron DOS will show signatures and artifacts of contact bondings [15], [24]. 

In the weak coupling (CB) limit however, integer charge addition is favored, and transitions 
between states that differ by a single electron appear as spectral signatures [26]. At equilibrium, 
it is convenient to introduce the Ground State Spectral Function by defining the Green's function 
in the Lehmann representation [21]: 

(N,0\c a \N + l,j)(N+l,j\cl\N, 0) 
E + i0+ — {Ef +1 - Eg) 

(iV,0|4|iV-l,j)(iV-l,j|c a |iV,0) 
E + i0+ - {Eg - Ef' 1 ) 

i[G(E)-&(E)] (7) 



G a p(E) = 
+ 

A aP {E) = 
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where a, (3 correspond to the orbital index, which in our case are the sites of the benzene 
molecule, and \N,i) denotes the i th excited state of a charge configuration of N electrons. The 
poles of this spectral function represent various transition energies for addition (removal) of 
electrons from the neutral ground state: 

Nr TTi JV jpN-1 

eg a = Ef +1 -E» (8) 
whose spectral strengths are given by: 



(N,0\c a \N + l,j)(N + l,j\cl\N,0) 



_Nr 

' Oj,aP 



T, 



Na _ /AT nUt 



0j,a/3 



(N,0\c}\N-lJ)(N-l,j\c a \N,0) (9) 



The first (addition) term adds an electron to orbital j3, taking the system from an N electron 
ground state to the j th (N+l) electron excited state, and then removes it from orbital a, bringing 
it back to ground state. The second (removal) equation first removes an electron from a and 
then adds it to /3. 

One can re-write the expression in terms of diagonal terms only, replacing the recurring index 
a with a single index, in a more convenient form as: 

<(£) = E [< r J(E ~ 4 r ) + r f a 5(E - e»;)} (10) 

3 

The spectral function shown in Fig. 2d represents the removal and addition strength of various 
transitions at their energies given by Eq. |T0j Notice that there are numerous peaks in this 
spectrum calculated from the many-electron transitions, due to the possible transfer to various 
excited states of charged species shown in Fig. 2d. It is important to note that although each 
transition has a non-trivial spectral weight given by Eq. they satisfy an overall sum rule that 
amounts to the total electron number in the system. We will see in subsequent sections that these 
transitions involving excited states show up directly as transport signatures frequently observed 
in experiments. 

III. Non-Equilibrium 

This section is devoted to the various unique transport signatures in the weak coupling (CB) 
regime, many of which have experimental significance. We elaborate on how various excited 
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states get accessed as a result of contacts maintained at different potentials (non-equilibrium), 
and what the SCF theory completely misses in this regime. Throughout this paper we describe 
the electrodes (contacts) using corresponding electrochemical potentials fiL and (j,r and coupling 
strengths 7^ and 7^. 



A. Coulomb Blockade approach: Rate equation model 

Transport in the CB limit [25]-[27] is often modeled with a rate equation approach, in which 
the steady state addition and removal of electrons is described with a rate equation for the 
nonequilibrium probability Pf of each N electron many -body state \N, i) with total energy Ef. 
The master equation involves transition rates R(N,i)^(N±i,j) between states differing by a single 
electron, leading to a set of independent equations defined by the size of the Fock space [26] 
dR N 



dt 



^2 [R(.N,i)->(N±l,j)Pi f ~ R(N±l,j)^{N,i)Pj* ±l (11) 



J 



along with the normalization equation J2i,N Pi — 1- We define rate constants 

V Nr = ry \ T Nr \ 2 

ija !a \ 1 ija I 

r^ a = t w Na \ 2 (\2) 

where 7 a represents lead molecule broadening or coupling via the end atoms, described using 
Fermi's Golden rule. These constants represent the partial probability for the electron to be 
injected by the end atom into a given many-electron ground or excited state. The transition rates 
are now given by 



R(N,i)^(N-l,j) - ^ija 1 _ f( e ij T - Va) 

a=L,R 

R(N-l,j)^(N,i) = ^ij r af( e ij T ~ Va)- (13) 

a=L,R 

for the removal levels (N,i <-> N — 1, j), and replacing (r — > a, f — > 1 — /) for the addition 
levels (N, % <-» N + 1, j). yt, a are the contact electrochemical potentials, / is the corresponding 
Fermi function, with single particle removal and addition energies e^ r = Ef — E^~ l , and 
efj a = E^ +1 — Ef. Finally, the steady-state solution to Eq.(fTTb is used to get the left terminal 
current as 

1 = ± ^ H [ R (N,i)^(N±i,j) P i N ~ R (N±i,j)^(N,i) P j N±1 _ (I 4 ) 
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where states corresponding to a removal of electrons by the left electrode involve a negative 
sign. 

B. STM Limit: 

We briefly elucidate the relationship between spectral functions defined earlier and STM 
conduction spectra. Electrical conduction depends on the measurement geometry [24] and charg- 
ing determined by capacitive voltage-division ratio rj between leads in the Laplace solution as 
opposed to the resistive voltage-division ratio 7 = 7_r/tl which determines the extent to which 
the levels are filled or emptied by the leads. The source and drain potentials are then given by 
Hl = E f + r]V d and fi R = E F - (1 - r))V d . 

Consider a simple picture in which one contact is very weakly coupled (7 << 1, 77 = 0), 
equivalent to the molecule being in equilibrium with left contact /i L . In the 77 = limit the 
molecular energy levels are pinned to this contact implying that for a positive voltage [ir < Hl, 
and hl remains at the equilibrium position. This picture is analogous to STM shell tunneling 
experiments [30], in which the weakly coupled STM tip acts as a voltage probe, thereby 
generating the single particle spectrum, the molecule/dot held in equilibrium with the more 
strongly coupled contact, in this case, the substrate. 

• Ground State Spectral function: It is expected, with a more strongly coupled contact, that 
the right contact voltage probe, such as an STM tip can add or withdraw into or out of 
the dot at energies corresponding to addition or removal energies defined in Eq. [T7] The 
stronger coupling to the left contact ensures that an electron be added or removed as soon 
as the tip removes or adds an electron thus maintaining overall charge neutrality. In this 
case conductance spectrum proportional to the equilibrium spectral function is obtained as 
shown in Fig. 3b. 

• Excited Spectral Functions: In the previous case, the chemical potential of the left contact 
is fixed above the transition level e^ r but below e^ r , thus maintaining the molecule's charge 
neutrality in its ground state (i.e., \N, 0)), and hence only the Ground State spectral signature 
AqAE) is observed. However, in a general non-equilibrium scenario, access to excited states 
of the neutral and charged molecule becomes feasible and hence description in terms of 
spectral functions corresponding to addition/removal from the i th excited state of the neutral 
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Fig. 3. STM limit - mapping spectral signatures: a) Removal Spectral functions of ground state Aq^/ r (E) (continuous) and 
first excited state A^ a (E) reproduced in the STM spectra in b). The STM spectra can also show signatures of charge neutral 
excited states shown dotted in a), depending on the position of the equilibrium chemical potential (see text), c) Simple schematic 
depicting the interplay of A^ r L ^ R {E) and A^ L / R (E) resulting in satellite peaks SI and S2 in conduction spectra of lower right 
half, d) Conductance spectra reproducing features of both ground and excited state spectral functions A^ r L ^ R {E) and A^ a (E). 



molecule is required: 

= E WJ{E - 4 r ) + r^6(E - eg-)] , (15) 

3 

where a now corresponds to the two sites that are coupled to the left (L)/ right (R) contacts. 

Nr 
10 



For example, let the equilibrium chemical potential be situated at a position above e N '' 



— Eq~ 1 , shown dotted in Fig. 3a. Given, a positive bias (/il > /ir) the above transition 
is energetically feasible only if the ground state of the cation (\N — 1, 0)) is accessed, which 
occurs for a tip voltage corresponding to hr below e^ r . Once this transition is accessed, 
spectral function Af L , R involving the first excited state gets involved due to the initial 
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condiction [i L > efj, due to which the neutral excited state \N,1) can be accessed. This 
results in additional satellite peaks SI and S2 in Fig. 3d. A schematic of transitions that 
consititute the satellite peaks SI and S2 due to A^ a (E) is shown in Fig. 3c. In this figure, we 
have only shown the removal levels for brevity, and extension of the argument by including 
addition levels is trivial. In general, in the STM regime, one can write a simple expression to 
evaluate the conductance formula as a weighted average over various excited state spectral 
functions: 

('« 

We have thus shown a simple signature that indicates the access of excitations in the many-body 
spectrum of the neutral molecule. In general, one may view near-equilibrium conduction in the 
CB regime using single particle energy levels: 

e% = E» - E?' 1 (17) 

and their corresponding spectral weights. 

C. Break Junction limit 

The break junction limit is achieved by setting 7] = 0.5, 7 = 1, implying that both contacts are 
equally coupled to the molecular dot and half the applied voltage appears across the molecular 
levels which in our case transition energies efj. The many-body configuration of the molecule 
consists of its ground state \N, 0) and the first excited state \N,1) separated by a gap similar to the 
HOMO-LUMO gap A, followed by a set of closely spaced excitations denoted by \N,i),i > 1. 
The I-V characteristics in this limit show certain key signatures which result from how these 
excitations are accessed. 

The onset of conduction is established by the offset between the equilibrium Fermi energy Ep 
and the first accessible transition energy e^ r . The qualitative shape of the I-Vs depends on how 
the excitations are accessed. Recall that 77 = 0.5 implies that the molecular levels are displaced 
with respect to the contact electrochemical potentials by the applied voltage. If the excited states 
are not accessed simultaneously or prior to the threshold transition e^ r , as shown in Fig. 4b, 
the I-V has a brief staircase of plateaus before a quasi linear rise in current. This quasilinear 
current rise occurs due to a huge number of closely spaced transport channels that are triggered 
only when transitions involving an excitation appear within the bias window. However, the 
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quasilinear current can also appear prematurely without an intervening plateau, if a feasible 
transition to an excited state appears in the bias window at or before the threshold transition. 
This situation is shown in Fig.4b, where e^ r also appears at threshold, resulting in a quasilinear 
regime immediately following the onset. The two distinct I-Vs have been observed experimentally 
[10]— [13], [29] and depend merely on the position of the equilibrium electrochemical potential 
(Fermi energy) with respect to the transition energies. In the meanwhile, similar SCF based 
I-V characteristics show adiabatically smeared out currents whose onsets get postponed by the 
changing position of the equilibrium E F , as shown in Fig. 4c. The SCF potential from Eq. 5, 
determines how levels float with respect to their non-equilibrium occupation [1]. It is readily 
seen by comparing Fig.3a,b with Fig. 3c that any self consistent potential cannot change the 
qualitative features of the I-Vs in order to resemble the CB features. 

D. Connection to experiments: Fitting Data using Coulomb Blockade model 

We consider matching the I-V shapes using our CB model using consistent fitting parameters. 
The experiments conducted on conjugated phenylenes [11], [12] at low temperatures (T = 30A') 
suggest strong Coulomb Blockade effects. It is worth noting that an 'orthodox' theory simply 
involving junction resistors and capacitors would also manage to capture the zero-bias suppressed 
conductance, the subsequent sharp onset and the linear current rise; however, it would not capture 
the intervening plateaus, fine structures in the I-Vs, and the gateability of the current levels and 
their asymmetry features that arise due to discrete transitions in the molecular configuration space 
[10], [14]. In contrast with metallic islands, a molecular dot shows significant size quantization 
that leads to quantum corrections to the junction capacitance, and gets further modified at high 
bias to involve nonlinear corrections to it arising from partial densities of states filled separately 
by the two contacts. 

While our model explains salient features of a lot of Coulomb Blockade experiments [10], 
[14], [29] it is interesting to note that in some cases the same molecule showed CB behavior at 
low temperature and SCF behavior at higher temperatures [11], [12]. A possible explanation is 
that at low temperature the molecule could be frozen into a configuration where the plane of the 
middle ring is oriented perpendicular to the side rings, while room temperature structures sample 
other configurations and are rotated on average. This is supported by the fact that current levels 
at room temperature are an order of magnitude greater, which can be attributed to an increased 
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Fig. 4. CB transport under identical contact coupling, a) Schematic of CB conduction resulting in qualitatively different I-V 
characteristics. For 1) /j,f = Mo one observes I-V with a Coulomb staircase with a plateau followed by a quasi-linear rise. 2) 
Hf = A*i one observes the quasi linear I-V upon reaching threshold, (b) This occurs with the intersection of /Uz, = euf an d 
fJ-R = e ocf l me m me stability diagram. Stability diagram shown for N = 6 particle blockade region, c) Distinct I-V's under 
cross sections [if — fio and /j,f — /ui. 

average degree of conjugation along the molecular backbone. In contrast at low temperature the 
rotated central ring has a weaker coupling with the rest of the backbone, which could reduce its 
broadening while increasing electron localization and charging, leading to CB behavior. In fact 
some of the experiments feature bulky middle groups like antracence. Steric side groups that are 
deliberately inserted to facilitate this rotation of the central rings and enforce CB [12]. While 
doing exact calculations on these molecular structures is beyond the scope of the present paper, 
we consider making simple fits by considering the following facts: 

• Current Levels: Using the fitting parameters 71 = 72 ~ 5 — lOmeV we obtain current levels 
similar to experimental data. It is important to note that changing 7 does not affect the 



February 6, 2008 



DRAFT 



16 





Fig. 5. Experimental fits for data [11], [12]. a) T = 30A", 7 « 5meV. b)T= 295A", 7 « 5meV c)T = 30K, 7 « 0.25meV. 



conductance before the threshold voltage which shows a vanishingly small pre-threshold 
current. 

• Threshold Voltage: We noticed in the last section with that the gap A between ground and 
first excited states of the neutral molecule is important in determining the qualitative shape 
of the I-V. When the equilibrium electrochemical potential E F lies above mid-gap between 
e^J' and e^ r , the first excited state becomes voltage-accessible before the ground state of the 
charged species is accessed and populated simultaneously via /j,r = €qJ, giving rise to the 
quasi linear I-V immediately following the very first current onset. In all the experimental 
data, we observe a threshold voltage between 0.5 — 0.7 V thus tuning the gap A w 0.6 — 0.8 
V. 
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Fig. Eta) (b) and (c) are fits obtained for experiments [11], [12]. In case of molecular asymmetries 
[11], only positive bias is considered [9], the I-V asymmetries themselves being attributed to 
polarization effects [7]. In obtaining an experimental fit for [12], in Fig.|5tb), we used T = 295 
K consistent with experiment. Notice that the first peak has broadened significantly. The higher 
temperature Coulomb Blockade could possibly be attibuted to the fact that molecule involved 
has an anthracene based middle ring that is much bulkier, thus leading to a higher temperature 
frozen configuration stabilized by steric interactions. 

The molecular system we consider is a simple prototypical molecule (benzene in our case) with 
calculations based on simple parameters that are associated with this minimal system. Performing 
calculations on a real molecule-electrode system will be needed to yield a quantitative fit in terms 
of threshold voltage, current levels and positions of peaks. However, the conduction mechanism 
remains the same. The exponentially larger configuration space of even a minimal Coulomb 
Blockaded molecule makes a first-principles calculation of its transport properties inordinately 
challenging compared to SCF treatments in the literature. However, the SCF calculations do 
not capture the non-equilibrium transition rates between the many-body states, which as we 
argued earlier carry crucial correlation signatures that are experimentally observable for ultrashort 
molecules. Such a "real" calculation involving the quantum chemistry of larger molecules and 
contact bondings within this nonequilibrium full CI treatment is still at a very early stage [31]. 
Furthermore, it needs to be supplemented with the broadening of the many-particle states that 
could affect the interference between nearby levels, an issue that has received relatively little 
attention [32]-[36] and requires further work. 

IV. CONCLUSION 

In this paper, we have developed a Coulomb Blockade approach for molecular conduction 
through short molecules using Benzene as a prototype. We have shown how equilibrium and 
non-equilibrium signatures are very different from the traditional NEGF-SCF viewpoint, and 
that the CB approach is appropriate in the weak coupling limit. Many I-V features distinct 
to the CB regime are often seen in experiments. These features that are easily obtained using 
a full Configuration Interaction master equation approach are potentially very hard to obtain 
within any effective one-electron potential, even for a minimal model. A particular challenge 
therefore lies in bridging the SCF and CB regimes while paying close attention to coherent 
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level broadening and associated interferences. The emergence of many recent experiments on 
molecular dots, exploring the interplay between charging, quantization and level-broadening, 
should prove invaluable in further theoretical developments in this regard. 
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